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Memorie deiia 




Abstract. We provide a summary of results, obtained from a multiwavelength (TeV y-ray. 
X-ray, UV, optical, and radio) campaign of observations of AE Aqr conducted in 2005 
August 28-September 2, on the nature and correlation of the flux variations in the various 
wavebands, the white dwarf spin evolution, the properties of the X-ray emission region, and 
the very low upper limits on the TeV y-ray flux. 
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1. Introduction 

AE Aqr is a bright {V ~ 11) novalike cat- 
aclysmic variable consisting of a magnetic 
white dwarf primary and a K4-5 V secondary 
with a long 9.88 hr orbital period and the short- 
est known 33 s white dwarf rotation period. 
Although originally classified and interpreted 
as a disk-accreting DQ Her star, AE Aqr is now 
widely believed to be a former supersoft X-ray 
binary (Schenker et al. 2002) and current mag- 
netic propeller (Wynn, King, & Horne 1997), 
with most of the mass lost by the secondary 
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being flung out of the binary by the magnetic 
field of the rapidly rotating white dwarf. 

Because of its unique properties and vari- 
able emission across the electromagnetic spec- 
trum, AE Aqr has been the subject of nu- 
merous studies, including the campaign of 
multiwavelength observations in 1993 October 
(Casares et al. 1996, and the series of pa- 
pers in Buckley & Warner 1995). Given 
the many improvements in observing ca- 
pabilities since that time, we undertook a 
campaign of multiwavelength observations 
of AE Aqr in 2005 August 28-September 
2, built around a Chandra/HST/YLA joint 
observing proposal. To these coordinated 
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Fig. 1. Top: Location of facilities participating in tiie 2005 multiwavelength campaign of observations of 
AE Aqr; stars mark the locations of CBA and AAVSO astronomers. Bottom: Timeline of the observations, 
with the name of the lead individual for each facility noted in parentheses. 



X-ray, UV, and radio observations, we added 
H.E.S.S. and MAGIC TeV y-ray observa- 
tions, Skinakas Observatory 1.3-m telescope 
and Mount Laguna Observatory 40-in tele- 
scope fast B-band optical photometry, and ex- 
tensive B- and V-band optical photometry ob- 
tained by Center for Backyard Astrophysics 
(CBA) and American Association of Variable 
Star Observers (AAVSO) professional amateur 
astronomers B. Allen, M. Bonnardeau (BZU), 
P. de Ponthiere (DPP), A. Gilmore (GAM), 
K. A. Graham (GKA), R. A. James (JM), 



M. Koppelman (KMP), B. Monard (MLF), A. 
Oksanen (OAR), T. Richards (RIX), D. Starkey 
(SDB), and T. Vanmunster (VMT) (see Fig. 1). 

As is maddeningly typical of multiwave- 
length campaigns, our efforts met with mixed 
success. H.E.S.S. was able to observe AE Aqr 
for only two brief intervals; the weather was 
less than ideal during the MAGIC observa- 
tions; our plan to obtain HST STIS time- 
tagged UV echelle grating spectroscopy was 
thwarted by the failure of the FUV-MAMA 
detector's power supply; our backup plan to 



Mauche et al.: Multi wavelength Campaign of Observations of AE Aqr 



-1 

I , , , , , I , , , , , , 
10 ^ Chandra: X — ray 



o = 



orbital phase 
1 




lO^'F 



3 10"^ r 



GALEX: NUV 



Ml \ 

li V 



10 - 

M 11 - 



Skinakas & Mt. Laguna: B — band Optical 



13 - 

I , , I , , , , I , , , , I , , , , 

^° 7 CBA & AAVSO: V- & B-band Optical 

™ 11 h ^^^^^-. ^v/iv^^ .^'U.-...^ 

13 - 



VLA: Radio 



6 10' r 




10" 



612.5 



613.0 613.5 
JD-2453000 



614.0 



Fig. 2. Chandra X-ray; GALEX NUV; Skinakas, Mount Laguna, CBA, and AAVSO B- and V-band (blue 
and red points, respectively) optical; and VLA radio light curves of AE Aqr. 
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Fig. 3. Chandra X-ray (upper histogram); GALEX NUV (blue points); Skinakas & Mount Laguna (black 
points) and CBA & AAVSO (red points) B-band optical; and VLA radio (lower histogram) light curves of 
AE Aqr. The times of some of the more prominent optical flares are marked with dotted blue vertical lines. 



4 



Mauche et al.: Multiwavelength Campaign of Observations of AE Aqr 



1000 1 20 




Chandra HETG Chandra HETG 

Fig. 4. Correlation between the Chandra HETG count rate and the GALEX NUV count rate {left) and 
the VLA flux (right). Linear fits to the data {blue lines) have correlation coefficients of +0.83 and -0.22, 
respectively. 



obtain GALEX FUV and NUV photometry was 
compromised by a then ongoing problem with 
the FUV detector; the timings of the pulsa- 
tions observed in the GALEX NUV data are 
uncertain because of an uncertain offset to the 
spacecraft clock; the observations were spread 
out over a longer interval and had less over- 
lap than is ideal, and some of the acquired 
data never have been analyzed. Nonetheless, 
we have published results to summarize and a 
number of new results to report. 



2. Light Curves 

The Chandra HETG (0.5-6 keV) X-ray, 
GALEX NUV (1750-2800 A) ultraviolet, 
Skinakas and Mount Laguna Observatory B- 
band optical, CBA and AAVSO B- and V-band 
optical, and VLA 3.6 cm radio light curves ob- 
tained during our campaign are shown in Figs. 
2 and 3. Although the flux is variable in ev- 
ery waveband, the nature of the variability in 
the radio waveband is clearly different from 
that in the optical through X-ray wavebands. 
The flares in these wavebands are highly cor- 
related, last between a few hundred and a few 
thousand seconds, and have an amplitude that 
increases from the optical through the NUV. 
To quantify the degree of correlation between 
the various wavebands, we show in Fig. 4 scat- 
ter plots of the GALEX NUV and VLA ra- 
dio flux as a function of the Chandra HETG 
X-ray flux. Linear fits to the data, shown by 
the blue lines, have correlation coefficients of 



+0.83 and -0.22, respectively, demonstrating 
a strong positive correlation between the X- 
ray and NUV wavebands and a weak negative 
correlation between the X-ray and radio wave- 
bands. We take these results as evidence that 
the optical through X-ray emission regions are 
tightly coupled and hence in close proximity, 
whereas the radio emission region is not. 



3. Spin Period Ephemeris 

de Jager et al. (1994) used optical and UV pho- 
tometry of AE Aqr, obtained over the inter- 
val 1978 June- 1992 November, to show that 
(1) the 33 s pulsation follows the white dwarf 
around the binary center of mass, with a pro- 
jected semi-amplitude {a sin/)/c = 2.04 + 0.13 
s, and (2) the rotation period of the white dwarf 
is spinning down at a rate P - 5.642(20) x 
10"''*. Mauche (2006) used the Chandra data 
obtained during our campaign to show that the 
33 s X-ray pulsation also follows the white 
dwarf around the binary center of mass, with a 
projected semi-amplitude (asinO/c = 2.17 + 
0.48 s. He also found that pulse timings de- 
rived from Chandra and from previous XMM- 
Newton (2001 November) and ASCA (1995 
October) observations are not consistent with 
the de Jager quadratic spin ephemeris, al- 
though the ephemeris could be fixed with the 
addition of a cubic term P = 3.46(56) x 
10"'^ d"'. Dulude (2009) derived pulse timings 
from three out of four nights of the Skinakas 
and Mount Laguna B-band optical photometry 
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Fig. 5. O - C residuals of fits to the optical timings of de Jager et al. (1994) [circles], the ASCA, XMM- 
Newton, and Chandra X-ray timings of Mauche (2006) [xs], and the Skinakas and Mount Laguna optical 
timings of Dulude (2009) [squares] for the de Jager et al. (1994) quadratic ephemeris (top) and the quadratic 
(middle) and cubic (bottom) ephemerides given in Table 1. Standard deviations of the data relative to the 
fits are shown by the error bars. Blue curve in the top panel is the additional P = 3.46 x 10"" d ' cubic 
term to the de Jager et al. (1994) ephemeris proposed by Mauche (2006). 



Table 1. Spin Ephemeris Constants 





To (BJD) 


P(d) 


P(d d"') 


P(d-') 


de Jager . 
Quadratic 
Cubic 


2445172.0000423(10) 
2445172.0000392(12) 
2445172.0000428(11) 


0.00038283263840(28) 
0.00038283263735(18) 
0.00038283263823(17) 


5.642(20) X 10-" 
5.752( 8) X 10-'* 
5.599( 7) X 10-"* 


4.97(31) X lO-" 



obtained during our campaign. He found that 
the optical pulse timings are consistent with the 
Chandra X-ray pulse time within the errors, 
and used all the available data to further inves- 
tigate the white dwarf spin evolution. We dis- 
covered a few errors in his transcription of the 
de Jager pulse timing data, so redid this anal- 
ysis for both quadratic and cubic ephemerides. 
The resulting fit residuals are shown in Fig. 5 
and the spin ephemeris constants are listed in 
Table 1, which includes the de Jager quadratic 
spin ephemeris constants for reference. 

In each case, the fits are unweighted and 
the error estimates are derived by assuming 
that the timing errors are given by the standard 
deviation of the data relative to the fit (hence. 



xl = 1). While the cubic ephemeris is not re- 
quired by the data, the spin-down rate P in 
the quadratic ephemeris differs from de Jager's 
value by S.lcr; regardless of how it is parame- 
terized, we confirm that the white dwarf in AE 
Aqr is spinning down at a rate that is faster than 
predicted by de Jager et al. (1994). 

4. X-ray 

We summarize here and in Fig. 6 the results 
obtained by Mauche (2009) from our Chandra 
HETG observation of AE Aqr First, the X- 
ray spectrum is that of an optically thin multi- 
temperature thermal plasma; the X-ray emis- 
sion lines are broad, with widths that increase 
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Fig. 6. Top: Chandra HETG spectrum of AE Aqr {black histogram), icr error vector {bhie histogram), 
and the best-fit model described in the text {red histogram). Middle: He-like //(/ + r) line ratios of O VII, 
Ne IX, Mg XI, and Si XIII. Blue rectangles delineate the Icr error envelope of the measured line ratio 
and inferred log for each ion. Bottom left: Spin-phase radial velocities of the X-ray emission lines {filled 
circles with error bars), best-fitting sine function {solid blue curve), and y velocity {dotted horizontal line). 
Bottom right: Orbit-phase delay in the arrival of the X-ray pulse maxima {filled circles with error bars), 
best-fitting cosine function {solid blue curve), and 2-s semi-amplitude cosine modulation observed in the 
optical {dotted cun-e). Adapted from Figs. 7, 8, and 10 of Mauche (2009) and Fig. 1 of Mauche (2006). 



with the line energy, from 510 km s ' for O 
VIII to 820 km s ' for Si XIV; the X-ray 
spectrum is well fit by a plasma model with 
a Gaussian emission measure distribution that 
peaks at logr(K) = 7.16, has a width cr - 
0.48, an Fe abundance equal to 0.44 times so- 
lar, and other metal (primarily Ne, Mg, and 
Si) abundances equal to 0.76 times solar; and 
for a distance c/ = 100 pc, the total emission 
measure EM = 8 x lO^"* cm""* and the 0.5- 
10 keV luminosity Lx = 1 x 10"" erg s"'. 
Second, based on the //(/ -i- r) flux ratios of 
the forbidden (/), intercombination (0, and re- 
combination (r) lines of the He a triplets of N 
VI, O VII, and Ne IX measured by Itoh et al. 
(2006) in the XMM-Newton RGS spectrum and 
those of O VII, Ne IX, Mg XI, and Si XIII 



in the Chandra HETG spectrum, either the 
plasma electron density increases with temper- 
ature by over three orders of magnitude, from 
He ~ 6 X 10"' cm-3 for N VI [log r(K) « 6] to 
He ~ 1 X lO'-* cm-3 for Si XIII [log r(K) ^ 7], 
and/or the plasma is significantly aff'ected by 
photoexcitation. Third, the radial velocity of 
the X-ray emission lines varies on the white 
dwarf spin phase, with two oscillations per spin 
cycle and an amplitude K ^ \6Q km s ' . 

5. TeV y-ray 

We summarize here the results obtained by 
Sidro et al. (2008) during four consecutive 
nights and a total of 15.3 hr of observations 
of AE Aqr with the 17-m MAGIC atmo- 
spheric Cherenkov telescope. Unfortunately, 
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on all four nights of observations, the cosmic 
trigger rate was low and the optical extinction 
was high, indicating that the sky was veiled by 
calima (airborne dust and sand, blown off the 
Sahara Desert). All of the MAGIC data were 
analyzed using the standard reconstruction and 
analysis software, and a set of special Monte 
Carlo simulations were employed to correct the 
energy scale and gamma efficiency of our anal- 
ysis. Because of the relatively poor observing 
conditions, the lack of a detection, and the ab- 
sence an entirely independent analysis of the 
data — because the analysis has not been thor- 
oughly vetted by the MAGIC consortium — 
the following results should be considered pre- 
liminary. Nevertheless: 

No evidence was found for steady y-ray 
emission on any night or during the sum of the 
four nights of our observations. On three of the 
four nights, the upper limits for steady emis- 
sion are ~ 8.0 x 10"'^ photons cm"^ s"' above 
340 GeV and ~ 2.5 x lO'^^ photons cm'^ s'^ 
above 650 GeV. Note that, using a total of 
68.7 hr of exposure of AE Aqr obtained with 
the Whipple Observatory 10-m atmospheric 
Cherenkov telescope, Lang et al. (1998) de- 
rived an upper limit for steady y-ray enaission 
of 4 X 10-1^ photons cm-^ s'^ above 900 GeV. 
These upper limits are three orders of magni- 
tude lower than the flux reported by others dur- 
ing the purported y-ray burst activity. 

In addition to steady emission, we searched 
for pulsed y-ray emission on each night and 
during the sum of the four nights of our ob- 
servations, using both the Rayleigh test and 
phase-folding on frequencies near the white 
dwarf spin frequency and its first harmonic. No 
evidence above 2cr was found for pulsed y-ray 
emission on any of the four nights, with an up- 
per limit of ~ 1.5 X 10"'^ photons cm"^ s"^ 
above 650 GeV. 

6. Summary and Conclusions 

• We interpret the strong correlation in the 
flux variations between the X-ray and NUV 
wavebands in AE Aqr as evidence that the 
flaring component of the optical through 
X-ray emission regions are tightly coupled 
and hence in close proximity. 



• We interpret the weak negative correlation 
in the flux variations between the X-ray 
and radio wavebands in AE Aqr as ev- 
idence that the flaring component of the 
radio emission region is not tightly cou- 
pled to and hence distinct from the optical 
through X-ray emission regions. 

• The close association of the optical through 
X-ray emission regions is demonstrated by 
the common phase and 2-s semi-amplitude 
of the binary-phase dependence of the 
pulse timings. More specifically, the pul- 
sating component of the optical through 
X-ray emission regions follow the white 
dwarf around the binary center in mass. 

• Although the fit parameters of the spin 
ephemeris of the white dwarf depend on 
how the pulse timings are modeled and 
how the pulse timing errors are handled, we 
confirm that the white dwarf in AE Aqr is 
spinning down at a rate that is faster than 
predicted by the quadratic ephemeris of de 
Jager et al. (1994). 

• The results of the Chandra HETG obser- 
vation of AE Aqr (Mauche 2009), specif- 
ically the high inferred plasma densities 
and the spin-phase radial velocity varia- 
tions of the X-ray emission lines, are in- 
consistent with the recent models of Itoh 
et al. (2006), Ikhsanov (2006), and Venter 
& Meintjes (2007) of an extended, low- 
density source of X-rays in AE Aqr, but in- 
stead support earlier models in which the 
dominant source of X-rays is of high den- 
sity and in close proximity to the white 
dwarf. 

• Given the non-detections of AE Aqr by the 
Whipple Observatory (Lang et al. 1998) 
and MAGIC (Sidro et al. 2008) atmo- 
spheric Cherenkov telescopes, and absent 
any publication of a positive TeV y-ray de- 
tection since the mid-1990's, we suggest 
the possibihty that, contrary to previous 
claims and common belief, AE Aqr is not 
a TeV y-ray source. 
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DISCUSSION 

ROBERT CONAN SMITH: Given the mea- 
sured radial velocity amplitude of the X-ray 
emission lines and the known spin period, can 
you estimate the radius at which the X-ray 
plasma is trapped? 



CHRIS MAUCHE: Evidently not. The mea- 
sured K velocity of the X-ray emission Unes 

is approximately 160 km s"', whereas the ro- 
tation velocity of material trapped on, and ro- 
tating with, the white dwarf is InrlP^^ia ~ 
1300(r//?wd) km s"^ A simple two-spot model 
of the X-ray emission region produces K as 
500 km s"^ , but this is still too large by a factor 
of ~ 3. This remains a puzzle. 

MANABU ISHIDA: You said that the esti- 
mate of the size of the X-ray emission region 
will be reduced by an order of magnitude rel- 
ative to that given by Itoh et al. (2006); isn't 
it still larger than the size of the white dwarf? 
Also, the maximum plasma temperature is too 
low for accretion onto the white dwarf. 

CHRIS MAUCHE: I actually said orders of 
magnitude: At the peak of the emission mea- 
sure distribution, the emission measure EM ^ 
5.8 X 10^^ cm"^ the density He ~ 1 X 10^"* cm"^ 
and the linear scale / = (EM/niy^^ w 1.4 x 
10*^ cm ~ 0.2 /?„d- I interpret the observed 
low maximum plasma temperature in AE Aqr 
as evidence that the accreting plasma does not 
pass through a strong hydrodynamic shock on 
its way to the white dwarf surface. The absence 
of such a shock might be expected for the low 
accretion rate, hence the low specific accretion 
rate, in AE Aqr. 

PIETER MEITJES: I disagree with your 
claim that AE Aqr is not a TeV source. We 
found in early TeV observations that the duty 
cycle is very small (<10%). So, many more 
observations, especially with H.E.S.S., are re- 
quired. 

CHRIS MAUCHE: Although I agree that the 
question of AE Aqr's purported TeV emission 
can be settled only with additional observa- 
tions, Whipple and MAGIC already have ob- 
served AE Aqr for a combined 84 hr and seen 
nothing. Either we have been extraordinarily 
unlucky, or the original claims that AE Aqr is 
a TeV source are wrong. I feel lucky. 



